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Abstract 


Petri nets are a general and well-established model of concurrent and 
distributed computation and behaviour, including that taking place in 
biological systems. In this survey paper, we are concerned with intrin- 
sic relationships between Petri nets and two formal models inspired by 
aspects of the functioning of the living cell: membrane systems and 
reaction systems. In particular, we are interested in the benefits that 
can result from establishing strong semantical links between Petri nets 
and membrane systems and reaction systems. We first discuss Petri 
nets with localities reflecting the compartmentalisation modelled in 
membrane systems. Then special attention is given to set-nets, a new 
Petri net model for reaction systems and their qualitative approach to 
the investigation of the processes carried out by biochemical reactions 
taking place in the living cell. 


Keywords: Petri net, biomodelling, membrane system, reaction sys- 
tem, locality, GALS, qualitative modelling, set-net, set membrane sys- 
tem. 


1 Introduction 


Petri nets (see e.g., [39]) are a general formal model for concurrent and dis- 
tributed computation. Over the years, an impressive variety of Petri net 
models suited for many different applications have been developed together 
with supporting theories and tools. A relatively new and ever more impor- 
tant field of application is biology. Thanks to its distributed character, the 
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Petri net approach appears to be particularly well-suited to provide com- 
putational and operational foundations for problems and issues arising in 
biology; see for example, [28], for a recent comprehensive overview of appli- 
cations of Petri nets in systems biology. On the other hand, to understand or 
make use of specific aspects of biological processes, new formal models have 
been proposed. Membrane systems and reaction systems are two examples 
of such models which are both abstractions of the functioning of the living 
cell. 


In this paper, we are concerned with intrinsic similarities and differences 
between Petri nets on the one hand, and membrane systems and reaction 
systems on the other hand. In particular, we are interested in establishing 
strong semantical links between these two models and Petri nets, and the 
possible mutual benefits that may result. Different enhancements of the 
Petri net model are considered for the faithful modelling of the dynamics 
of the biological phenomena represented by membrane systems and reaction 
systems. It is our aim to demonstrate the fruitful two-way interaction be- 
tween Petri nets and the other two models. We will be interested in Petri 
net semantics which open the way to importing Petri net methodologies and 
tools to the two biologically motivated computational models. Such seman- 
tics should be faithful, so rather than giving a Petri net interpretation, we 
adapt and incorporate new concepts into the Petri net framework, while re- 
taining the underlying Petri net philosophy. In particular, we will discuss 
an extension of the standard PT-nets with a concept of locality, and much 
emphasis will be given to the completely new class of Petri nets called SET- 
nets. The latter are suited for qualitative rather than quantitative modelling 
which plays an important role in rendering of the biochemical processes that 
take place in living cells. 


Petri nets are a graphical modelling language with strong mathematical, 
algorithmic and tool support for the specification and analysis of distributed 
systems. Many different classes of Petri nets have been developed since their 
first appearance in [37]. Their main common underlying philosophy is that 
states are distributed and actions have a local cause and effect (on the ad- 
jacent components of the net); for more discussions on this see |7]. The 
most typical Petri nets are without doubt the Place/Transition nets (or PT- 
nets) [8]. They are based on the production and consumption of resources 
by actions taking place in the system. Hence the resulting computational 
processes are essentially multiset based. Another well known, more funda- 
mental, Petri net class are the Elementary Net systems (or EN-systems) [8]. 
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Their dynamics is based on holding or not holding of local conditions rather 
than being resource based. 


Like PT-nets, membrane systems ([33, 34, 35, 36]) are essentially mul- 
tiset rewriting systems. As a computational model they are inspired by the 
way chemical reactions take place in cells which are divided by membranes 
into compartments. The reactions are abstracted to rules that specify which 
and how many molecules can be produced from given molecules of a certain 
kind and quantity. The dynamic aspects of the membrane system model 
including potential behaviour (computations), derive from such evolution 
rules. To capture the compartmentalisation of membrane systems, PT-nets 
are extended with transition localities. This makes it possible to have lo- 
cally synchronised executions, but it requires an extension of the causality 
semantics of PT-nets. 


Reaction systems [3, 10, 11, 12] are also a model for the investigation of 
processes carried out by biochemical reactions in living cells. The model is 
meant to contribute to the understanding of the interactions between such 
reactions. This time, however, biochemical reactions are based on qualitative 
rather than quantitative presence of resources. Hence, in order to obtain 
a faithful Petri net representation of reaction systems, it is necessary to re- 
evaluate the existing Petri net modelling approaches. We therefore introduce 
a new class of Petri nets, called SET-nets, that supports set-based (boolean) 
operations on tokens rather than the standard Petri net multiset-based token 
manipulation. 


Finally, we bring together the qualitative approach of reaction sys- 
tems and the compartmentalisation of membrane systems as the multiset 
approach of membrane systems is not always realistic from the point of view 
of explicitly counting huge number of molecules and reactions. Moreover, 
the resulting infinite state space makes it impractical or impossible to apply 
formal verification techniques. Therefore, we propose to consider set mem- 
brane systems, that is membrane systems with qualitative evolution rules. 
This is especially attractive as localities and SET-nets can be combined to 
yield a satisfactory Petri net semantics for set membrane systems. 


In this survey paper, we mainly describe approaches and give the essence 
of key results. More details on Petri nets and membrane systems can be 
found in [23]. Set membrane systems were introduced in [24], and SET-nets 
were first presented at the BioPPN Workshop held in Newcastle upon Tyne 
in June 2011 (see |27| for the informal workshop version). 
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2 Preliminaries 


Multisets A multiset over a set X is a function py: X > N = {0,1,2,...}. 
(In this paper we only consider the case that X is finite.) Multiset py is said 
to be empty if there are no x such that x € yw by which we mean that x € X 
and p(x) > 1. The empty multiset is denoted by 2. 

A multiset may be represented by listing its elements with repetitions, 
e.g., W = {y,y, z} is such that u(y) = 2, uw(z) = 1, and p(x) = 0 otherwise. 
We treat sets as multisets without repetitions. 

For two multisets « and py’ over X, the sum 4+ py’ is the multiset given 
by (wt p’)(x) = p(x) + p(x) for all xe € X, and if k € N then k- p is the 
multiset given by (k- w)(2) = k- p(x) for all « € X. The difference ps — pi 
is given by (wu — p’)(x) = max{p(x) — p’(x),0} for all ce € X. We denote 
pu < pw’ whenever p(a) < p’(x) for all ex € X, and uw < py’ whenever p< pi 
and spc Fp’ 


Petri Nets A Place/Transition net (or PT-net) is defined as a tuple 
PT = (Pl, Tr, W, Mo), 


where: Pl and Tr are finite disjoint sets of respectively places and transitions; 
W : (Tr x Pl) U(PIx Tr) > N is the arc weight function; and Mp: Pl > N 
is the initial marking (in general, any multiset of places is a marking). 
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Figure 1: PT-net. 


In diagrams, like that in Figure 1, places are drawn as circles, and 
transitions as boxes. If W(x, y) > 1, then (x,y) is an arc leading from x to 
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y. An arc is annotated with its weight if the latter is greater than one. A 
marking M is represented by drawing in each place p exactly M(p) tokens 
(small black dots). 

A step U of PT is a multiset of transitions. Its pre-multiset and post- 
multiset of places, °U and U°, are respectively given by 


*U(p) = J“ U(t)-W(p,t) and U*(p) = (Ut) -Wtp) , 
teU teU 


for each place p. For the PT-net in Figure 1 we have: 
epi 7Ti 731i, — 7b ~b _a rll ,T11 ,7311® _f,a@ ,a ,c ,a a 
{717,717,737} = {a1, 71, m3} and {777,717,737 }° = {at wT, T1, 13, 73 }- 


We can distinguish three modes of execution for PT-nets (from sequen- 

tial to fully synchronous). To start with, U is free-enabled at a marking M 
if °U < M. A free-enabled U is then: min-enabled if |U| = 1; and maz- 
enabled if U cannot be extended by a transition to yield a step which is 
free-enabled at M. For the PT-net in Figure 1 we have that, at the initial 
marking Mo, the step {7;1*, 7571} is free-enabled, {7;1!} is min-enabled, and 
{770,777,757 7°? } is max-enabled. 
That is, U is free-enabled at M if in each place there are enough tokens for 
the specified multiple occurrence of each of its transitions (note that each 
transition t needs to consume from each place p exactly W(p, t) tokens which 
cannot be shared with any other transition). Interleaving (min-enabledness) 
allows only one transition to be executed at a time. Maximal concurrency 
(max-enabledness) means that extending U would demand more tokens than 
M supplies. 

For each mode of execution m € {free, min, max}, a step U which is 
m-enabled at a marking M can be m-ezecuted leading to the marking M’ 
given by 

M'=M-—*U+U°. 
We denote this by M[U)mM’. Moreover, an m-step sequence is a finite or 
infinite sequence of m-executions starting from the initial marking. For the 
PT-net in Figure 1 we have 


Molrt a3} eed a mt oS 5a : 
Petri Nets with Inhibitor and Activator Arcs A PT-net can be 


equipped with two other kinds of arcs which test for the presence or ab- 
sence of tokens in places. More precisely, 


Inh C Pl x Tr and Act C PI x Tr 
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are respectively the sets of inhibitor and activator arcs. In diagrams, an 
inhibitor arc ends with a small open circle, while an activator arc ends with 
a small black circle. The role of both kinds of test arcs is to constrain the 
enabling of a step U by stipulating that it is free-enabled at a marking M if 
°U <M as well as 


e p¢ M whenever there is t € U such that (p,t) € Inh 
e p © M whenever there is t € U such that (p,t) € Act. 


All the remaining notions are the same as for PT-nets. 


Petri Nets with Localities A Place/Transition net with localities (or 
PTL-net) is defined as a_ tuple PTL = (Pl, Tr,W,£,Mo) such that 
(Pl, Tr, W, Mo) isa PT-net and @: Tr > Nisa locality mapping. In diagrams, 
such as that in Figure 4, boxes representing transitions belonging to the same 
localities are displayed on a grey background of the same shade. 

Localities can be used to define one more kind of enabling for steps of 
transitions. A step U of PTL is lmaz-enabled if U cannot be extended by 
any transition t satisfying @(t) € @(U) to yield a step which is free-enabled 
at M. That is, locally maximal concurrency (Imax-enabledness) is similar to 
maximal concurrency, but now only active localities cannot execute further 
transitions. For the PTL-net in Figure 4 we have that {7;'',7;'*} is Imax- 
enabled at the initial marking, but {7/"'} is not. All the remaining notions 
are the same as for PT-nets. 


3 Membrane Systems and Petri Nets 


A membrane structure ys (of degree m > 1) is given by a rooted tree with 
m nodes identified with the integers 1,...,m. We will write (7,7) € mw or 
i = parent(j) to mean that there is an edge from i (parent) to j (child) 
in the tree of 4, and i € yz to mean that 2 is a node of pw. The nodes 
of a membrane structure represent nested membranes which in turn deter- 
mine compartments (compartment j is enclosed by membrane j and lies 
in-between j and its children, if any), as shown in Figure 2. 

Let V be a finite alphabet of names of objects (molecules). A basic 
membrane system over jp is a tuple 


T= (V,p,w),-.-, we, Ri, <.-; Rm) 
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Figure 2: A membrane structure (m = 3) and its compartments with 1 
being the root node, (1,2) € w and 1 = parent(3). 


such that, for every membrane i, w? is a multiset of objects, and R; is a 
finite set of evolution rules r of the form ths’ > rhs", where lhs’ (the left 
hand side of r) is a non-empty multiset over V, and rhs” (the right hand 
side of r) is a non-empty multiset over 


VU {Gout |@eV}U {ain, |a € V and (i,j) € wu}. 


Note that a symbol ain; represents an object a that is sent to a child node 
(compartment) 7 and ao, means that a is sent to the parent node. If 7 is 
the root of yw then no indexed object of the form ao, belongs to rhs”. A 
configuration of II is a tuple 


C = (wy,.--;Wm) 


of multisets of objects, and Co = (w?,...,w?,) is the initial configuration. 
Figure 3 shows a basic membrane system over the membrane structure from 
Figure 2. 


A membrane system evolves from configuration to configuration as a 
consequence of the application (or execution) of evolution rules. There is 
more than one strategy in which this can be done. Maximal concurrency used 
to be the standard execution mode for membrane systems. Later, however, 
also in view of the intrinsic connections with Petri nets, other execution 
modes attracted interest. In particular, with the concept of localities added 
to PT-nets to represent compartments, locally maximal concurrency came to 
light as a new realistic execution semantics for membrane systems. Hence, 
similarly as in the case of PTL-nets, we can distinguish four such execution 
modes, all based on the notion of a vector multi-rule. 

A vector multi-rule of II is a tuple 


r= (11,...,lm) 
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Figure 3: A basic membrane system Io. 


where, for each membrane 7 of fu, r; is a multiset of rules from R;. For such 
a vector multi-rule, we denote by lhs; the multiset 


Ss" r;(r) - dhs” 


reER; 


in which all objects in the left hand sides of the rules in r; are accumulated, 


and by rhs; the multiset 
S> ri(r) + rhs" 


reR; 


of all (indexed) objects in the right hand sides. The first multiset specifies 
how many objects are needed in each compartment for the simultaneous 
execution of all the instances evolution rules in r. 

A vector multi-rule r of II is free-enabled at a configuration C if lhs; < 
w;, for each i. A free-enabled r is: min-enabled if |r1] +--+ + |rm| = 1; 
maz-enabled if no r; can be extended to yield a vector multi-rule which is 
free-enabled at C; and Imaz-enabled if no non-empty r; can be extended 
to yield a vector multi-rule which is free-enabled at C’. For example, in 
Figure 3, 


e (2,2, {r31}) is not free-enabled; 
e ({rii},@,@) is min-enabled but not Imax-enabled; 


e ({r11, 712}, 2, @) is lmax-enabled but not max-enabled; and 
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e ({ri1, 712}, {r21, 722}, @) is max-enabled. 


If r is free-enabled (free) at a configuration C, then C has in each 
membrane 7 enough copies of objects for the application of the multiset of 
evolution rules r;. Maximal concurrency (maz) requires that adding any 
extra rule makes r demand more objects than C’ can provide. Locally maxi- 
mal concurrency (/maz) is similar but in this case only those compartments 
which have rules in r cannot enable even more rules; in other words, each 
compartment either uses no rule, or uses a maximal multiset of rules. Mini- 
mal enabling (min) allows only a single copy of just one rule to be applied 
any time. 

The effect of the application of the rules is independent of the mode 
of execution m € {free, min, maz, lmax}. A vector multi-rule r which is m- 
enabled at C can m-evolve to a configuration C’ = (wj,...wy/,) such that, 
for each i and object a: 


w;(a) = w;(a) = Ihs; (a) oa rhs; (a) zo Ths arent(i) (Gini ) a > rhs} (Gout) 
i=parent(j) 


where Ths yyarent(i) = @ if i is the root of yz. We denote this by C >, C’. 
Moreover, an m-computation is a sequence of m-evolutions starting from the 


initial configuration. For the example in Figure 3 we have: 
Co in ({a, b}, {a, b, Cc, C, cH, {a}) ’ 


where r = ({r11,r12}, S, @). 


3.1 Petri Net Modelling of Membrane Systems 


There is a natural way of translating a basic membrane system II into a be- 
haviourally equivalent PTL-net PTL(II) = (P,T,W, 2, Mo), where multisets 
of places are used to represent the availability of molecules within the com- 
partments, and transitions correspond to evolution rules. Each transition is 
associated with a compartment and this information is represented by the 
localities of net transitions. The constructed PTL-net PTL(II) has a separate 
place mS for each molecule a and membrane j, and a separate transition 7)" 
with locality i for each rule r in compartment 7. The initial marking inserts 
w;(a) tokens into each place nm. The connectivity between transition t = 7; 
and place p = Ts is given by: 
if ghee 38 
woes { ce (ay ih 7 


otherwise . 
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as well as: 

rhs" (a) ey 

rhs" (Gout) if 7 = parent(i) 
rhs"(ain;) if i = parent(j) 
0 otherwise . 


W(t, p) = 


Figure 4 shows the result of the above translation for the basic membrane 
system in Figure 3. 
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Figure 4: PTL-net PTL(IIp) modelling the basic membrane system IIo. 


The pTL-net PTL(II) provides a faithful representation of the behaviour 
of the basic membrane system II. To capture this very close relationship, we 
define two bijective mappings, v and p, which allow us to move between II 


and PTL(I): 


e For every marking M of PTL(II), v(M) = (wi,...,Wm) is the config- 
uration of II, given by wij(a) = M(z%), for every molecule a. 


e For every step U of PTL(ID), p(U) = (r1,...,%m) is the vector multi- 
rule of II, given by r;(r) = U(z77), for every compartment 7 and every 
rule r € Rj. 


It is then possible to establish a direct relationship between (the operation 
of) the original membrane system and the PTL-net resulting from the above 
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translation at the system level: 


Com OT => vX(C) [p71 @))m YC!) (1) 
1 
M[U)mM’ => v(m)“, v(t’) 


for all modes of execution m € {free, min, maz, lmax}, configurations C of II 
and markings M of PTL(II). Together with v(Mo) = Co, such a result means 
that the (finite and infinite) m-step sequences of PTL (II) faithfully represent 
m-computations of II, and the same applies to markings and configurations. 


3.2 Petri Net Analysis of Membrane Systems 


Thanks to the very tight behavioural correspondence between II and PTL(II) 
captured by (1) above, analytical techniques developed for Petri nets can be 
applied to membrane systems. For example, one can use the invariant analy- 
sis based on linear algebra [42] to verify properties of configurations reachable 
from the initial one. If we take the membrane system in Figure 3 and apply 
the invariant analysis to the corresponding Petri net in Figure 4, then one 
can deduce that the total number of molecules a and b in compartment 2 is 
constant, irrespective of the initial configuration. 

Another direction is to use the causality semantics approach of Petri 
nets based on occurrence nets, allowing one to analyse entire computations 
rather than individual reachable configurations. In particular, occurrence 
nets allow one to investigate causality, concurrency and executability in 
system behaviour. In case of Imax-step sequences, however, one needs to 
modify the standard process construction, as first outlined in [26]. 

To analyse the state space of II one can also employ the reachability 
graph of PTL(II). Investigating reachability graphs has a long tradition in 
the field of Petri nets, and has produced several fundamental results. For 
example, reachability for PT-nets is decidable |[31, 29] which means that the 
problem of deciding whether a basic membrane system has a free- or min- 
execution leading to a given configuration can be decided, even when there 
are infinitely many reachable configurations. Another relevant property of 
II is whether the concentration of specific molecule(s) in specific compart- 
ment(s) can grow unboundedly. This problem, known in the are of Petri nets 
as boundedness, can be tackled using the coverability tree construction [18]. 
Coverability trees can also be used to decide whether two specific molecules 
can ever be simultaneously present in the same compartment. 
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4 Reaction Systems and Petri Nets 


Reaction systems [10, 11, 12] are a formal framework for the investigation 
of processes carried out by biochemical reactions. Thus the framework is 
inspired by biochemistry and its underlying ideas are motivated by the facil- 
itation/acceleration and inhibition/retardation, properties shared by a great 
number of biochemical reactions. Reaction systems constitute a computa- 
tional approach inspired by nature and are targeted at the investigation of 
ongoing dynamic changes occurring in biochemical systems through infor- 
mation processing. However, the model is based on principles remarkably 
different from those underlying other existing models of computation. 
A reaction system is a pair 


A= (S,A), 


where S is a finite background set comprising the entities of A, and A is the 
set of reactions of A. Each reaction is a triple of the form a = (R,I, P), 
where the three components are finite sets:* 


e RCS is the set of reactants; 
e I CS is the set of inhibitors; and 
e PCS is the set of products. 


The components of a reaction a = (R,J,P) may be denoted, respectively, 
by Ra, Ig and P,. 

A state of a reaction system is any set C' of its entities. Then an 
initialised reaction system is a triple 


A= (S,A,Co) , 


where (5, A) is a reaction system and Co C S is the initial state. A reaction 
system with background set S has exactly 2!5! potential states. To describe 
possible moves between these states, we need to say what is meant by an 
occurrence of a reaction or a set of reactions. 

A reaction a is enabled at a stateC C Sif Rg CC and IgnC = @; 
moreover, for the purpose of establishing the relationship with Petri nets, 
in this paper we will say that a set R of reactions is enabled at C' if each 


4 In the original definition these sets are assumed to be non-empty and RNI = @. 
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reaction of R is enabled. In such a case, R can occur with its effect on C 
being given by 


resR(C) = U Pes 
aceR 


We denote the resulting state change by C ma’ resr(C). If R is the set 
of all reactions enabled at C, then we may simply write C —> C’, where 
C! = resp(C). 

In the state change as described above, all the entities in C \ Ujer Pa 
disappear when 7? occurs. As a result, and unlike in other formal models of 
dynamic systems, there is no persistency in a reaction system in the sense 
that an entity present in a state disappears unless it is sustained by at least 
one reaction in . Consider, for instance, an initialised reaction system 


Ao = ({q, r, s}, {a, b, Ch {a}) ’ 


with background set {q,7r,s}, the initial state {q}, and three reactions: 


a = ({r, q}, 2, {r}) b= ({a}, {s}, tr, ah) c= ({g}, 2, {s}) - 


Then we have the following examples of state change: 
{a,c {5} {a} 
{r,q,5} =F frs} {rab irak Cars} Sh fry, 


One may observe that there is no conflict between reactions in the sense 
that the occurrence of one reaction might imply that another reaction which 
is also enabled at the current state, cannot occur. This, again, is a feature 
not found in most other formal models of dynamic systems. 

It is crucial to point explicitly to the ‘non-counting’ features of reac- 
tion systems: entities are either present or not, and produced or not, and 
reactions can or cannot occur given the presence or absence of certain enti- 
ties. There is no representation of multiple instances of entities or multiple 
occurrences of reactions. 

In general, reaction systems may have an environment and then operate 
within a changing context (with entities coming from the environment at 
each stage of evolution). Here, however, we will consider context-independent 
processes defined by a reaction system with an initial state provided by the 
environment, and every next state obtained as the result of reactions taking 
place in the previous state. 
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4.1 Set-Nets 


In [27| we investigate how to construct Petri net representations of reaction 
systems. While doing so, we made some general observations and assump- 
tions about the relationship between reaction systems and nets. 


e Entities can be represented by places, and reactions by net transitions. 


e Since there are no conflicts between reactions, activator arcs can be 
used to test for the presence of reactants (rather than claiming re- 
sources for the exclusive use as with ordinary arcs and input places). 


e Inhibitor arcs can be used to test for the absence of reactants. 


e All reactions that can occur in a reaction system do occur, and the 
only entities left after a state change are the newly generated products. 
In the Petri net framework, these features correspond to maximal con- 
currency and place resetting implemented by reset arcs [9]. 


We tried four different modelling methods, including high-level Petri nets |17]. 
In each case, we established a close correspondence between the evolutions of 
two corresponding models. All these net models, however, exhibited deficien- 
cies w.r.t. simplicity and/or elegance and/or tractability of the translation. 

In particular, in all four cases, one state of a reaction system would 
correspond to many markings of a corresponding Petri net, which is dramat- 
ically different from the one-to-one relationship between the configurations 
of a membrane systems and markings of a the corresponding PTL-net de- 
scribed in Section 3. We therefore proposed a new class of Petri nets, called 
SET-nets, which provide a stronger match with reaction systems and their 
semantics. 

The main idea is that in a SET-net there is no concept of counting. 
Places are marked or not marked and arcs have no weights. Set-nets resem- 
ble elementary net systems (EN-systems) [38] which is a fundamental model 
to study basic features of concurrent systems, including conflict, causality 
and independence. However, their execution semantics is different. In SET- 
nets, a marked place indicates the presence of a resource without any quan- 
tification. Hence any number of transitions that take input from this place 
can be fired at the same time. Moreover, firing a transition empties all its 
input places. Thus there are no conflicts over tokens in SET-nets, unlike in 
EN-systems or PT-nets. Similarly, places do not count the tokens, and the 
firing of a transition simply marks each of its output places (whether or not 
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they were already marked). We will build up the new model in two stages, 
introducing first SET-nets with only flow arcs. 
A SET-net is a tuple 


SN = (Pl, Tr, W, Mo) 


such that the four components are as in the definition of PT-nets, under the 
proviso that W always returns 0 or 1, and the initial marking Mp is a set (in 
general, markings are now sets of places). The firing rule for SN assumes 
that each step U is a set, and we denote by °U and U® the sets of all places 
p such that there is a transition t € U with W(p,t) = 1 and W(t,p) = 1, 
respectively. We then say that U is enabled at a marking M if °U C M. 
In such a case, U can be executed with its effect on M being given by the 
resulting marking 


M'=(M\*U)UU*. 


We denote this by M[U)M’'. Moreover, if U is the set of all transitions 
enabled at M (i.e., all transitions t such that {t} is enabled), then we ob- 
tain a maximally concurrent execution denoted by M[U)macM’' or simply 
M|) mar M’ as the maximally concurrent step enabled at M is unique. Hence 
a step U enabled at a marking M may contain two distinct transitions t and 
u for which *tN°u 4 @ or t?Mu® ¥ @ and yet the common places will never 
contain more than one token. 

To model reaction systems, we need additionally inhibitor arcs to cap- 
ture the effect of inhibitors in reaction rules. We therefore consider SET-nets 
with inhibitor arcs Inh C Pl x Tr. In such a case, the enabling relation 
changes similarly as for PT-nets with inhibitor arcs, and we say that a step 
U is enabled at a marking M if 


°UCM and (MxU)NInh=2@. 


The result of executing an enabled step remains the same as before. 


4.2 SEtT-Nets Modelling of Reaction Systems 


Reaction systems and SET-nets fit together well in the sense that both do not 
count tokens and both change states on the basis of the presence/absence 
of resources, represented by sets. Moreover, under the SET-net semantics, 
ordinary arcs (transitions) can be used to empty places. Finally, following 
the assumption that all reactions that can take place do take place, the 
maximal set-semantics can be employed. 
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Figure 5: SET-net SN(Ap) modelling reaction system Ag. 


Modelling inhibition aspects of reactions is rather straightforward, as 
illustrated by the SET-net in Figure 5, representing the context-independent 
initialised reaction system Ag considered earlier. As we already mentioned, 
places represent entities. Transitions 77, 7° and r° represent reactions, and 
r|, q and s, ensure that once the SET-net is active only tokens produced in 
the last maximal step are present in the current marking. Using inhibitor 
arcs gives a compact translation of reaction systems which is in a sense 
minimal w.r.t. the number of places, arcs and transitions. Moreover, relating 
the behaviour of the resulting SET-nets and the original reaction systems can 
be done as before. 

Formally, the places, transitions and initial marking of the translation 
SN(A) are given by: 


PE = {x*|seS} 
Tr = {r*|aeE A}UL{l,| s € S} 
Mo = {r*|sECo}. 


The flow and inhibitor arcs are as follows: 


W = {(n°,ls)|s eS} U 
{(m*°,7°) |aE AASE Ra} U 
{(7*,7°) |ae AAs eE P,} 

Inh = {(8*°,7%)|a@eEAAsEe ly}. 


Relating the behaviour of the SET-net model SN(.A) and the original 
reaction system A is straightforward using the mappings v(M) = {s | 7° € 
M}, for every marking M C Pl, and p(U) = {a | r* € U}, for every step 
U C Tr. 

It is then possible to establish a direct relationship between (the oper- 
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ation of) reaction systems and SET-nets at the system level: 


Cc4.c’! = WW: pU)=RAv-"C) [V) v1UC) 
. (2) 
Mu)M = v(M)Sum) 


for each state C of A, and each marking M of SN(A). Together with 
v(Co) = Mo, such a result means that the (finite and infinite) step sequences 
of SN(A) faithfully represent computations of A, and that there is a one-to- 
one correspondence between states and markings. 


Moreover, the maximally concurrent semantics of the SET-net corre- 
sponds to the execution of the reaction system, i.e., the execution of all 
enabled reactions in each evolution step: 


C—->C' = v-"(C) heere: v-\(C") 
MD moxM! —> v(M)—>v(M") 


for each state C of A, and each marking M of SN(A). 


Note that the fundamental class of EN-systems [38] extended with in- 
hibitor as well as activator arcs [16, 30, 32] basically has the same static 
structure as SET-nets. However, their treatment of conflicts between transi- 
tions accessing the same token, as well blocking a transition which could add 
a token to a marked place, are totally different. The latter issue has been 
noted in the past, and the constraint relaxed. For example, there are varia- 
tions of Petri nets, such as Boolean Petri nets, where adding a token to an 
already marked place does not add another token [5, 6, 14]. Also, behaviour 
of this kind was mentioned in |2] in the context of net synthesis. Having said 
that, the semantics considered in prior works was based on single transition 
firings, rather than (maximal) steps as is the case for SET-nets, and so the 
issue of ‘token sharing’ was never explicitly considered. 


The main initial motivation of our investigation was to see how Petri 
net based concepts could be deployed to analyse reaction systems. In doing 
so, we introduced the model of SET-nets which is an original contribution to 
the field of Petri nets. In the next section we will see how the introduction 
of SET-nets has motivated the introduction of a new model of membrane 
systems. 
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rir: {b} > {a} l 
{a, b} r12: {a} > {b, cing, Ging } 
ris? AOh —- 1ch 
2 3 
{a, b,c} i) 
ro1: {a,c} — {b} 731: {a} > {a, cour} 


ro2: {b} > {a} 


Figure 6: A basic set membrane system Xo. 


5 Petri Nets and Set Membrane Systems 


We now discuss membrane systems which use ‘qualitative’ rather than ‘quan- 
titative’ application of evolution rules to change the current state. The for- 
mal definitions and representation for this class of membrane are as those 
in Section 3, except that we are now working with sets rather than multi- 
sets of objects and evolution rules (similar to the operation of membrane 
systems in [1] where a qualitative approach was used in the application of 
rules within membranes, but a quantitative on for sending objects to the 
environment). 

A basic set membrane system over the membrane structure p is a tuple 


B= (Viet ce ot Risaiey Ri) 


as in Section 3, where each w? is a set of objects, and the left and right hand 
sides of every evolution rule are non-empty sets. Similarly, each configura- 
tion is composed of sets of objects. 

A vector set-rule of © is a tuple 


r= (11,...,lm) 


where, for each membrane 7 of j4, r; is a set of rules from R;. For two vector 
set-rules, r and r’, we denote r <r’ ifr; C r’;, for each i < m; andr <r’ if 
rCr’andr¥r’. For a vector set-rule r and i < m, we respectively denote 
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by: 
lhs; = U hs” and rhs; = U rhs" 


rer; rer; 


the set of all the objects in the left hand sides of the rules in r;, and the set 
of all the (indexed) objects in their right hand sides. We then say that a 
vector set-rule r is free-enabled / min-enabled / max-enabled / lmax-enabled 
exactly as in Section 3. Following this, a configuration C = (w1,...Wm) can 
m-evolve by a vector set-rule r which is m-enabled at C,, to a configuration 
C’ = (w},... wi) such that, for each compartment i of pu: 


m 


wi = WwW; ‘ lhs; U {a EV | ae rhs; V Qin; € rhs’, (2) m 


parent 


A(t, j) € ft Gout € rhsjz} . 


We denote this by C —>m C’. An m-computation is then defined as a (finite 
or infinite) sequence of consecutive m-evolutions starting from Co. 

The difference between the ‘qualitative’ and the ‘quantitative’ interpre- 
tation of the evolution rules is twofold. First, there may be two enabled 
evolution rules in a compartment with a common object in their left hand 
sides while there is only a single representant of that object in the current 
state in the compartment. In the current qualitative set-up, the two rules can 
be executed together. Second, if two simultaneously executed rules produce 
the same object in the same compartment, instead of adding two instances 
of this object, only one is added (so that we never have more than a single 
representant of an object in any given compartment). As a consequence, 
there is no need to use multisets of objects present in any single compart- 
ment to represent the current state, and there is no need to use vectors of 
multisets of rules in set membrane systems. In either case, using sets is fully 
sufficient. One may observe that with this view of state representation and 
system execution, maz-evolution is deterministic in set membrane systems. 


5.1 SEtT-Nets Modelling of Set Membrane Systems 


To faithfully capture the behaviour of basic set membrane systems, we need 
to extend SET-nets with localities. A SET-net with localities (or SETL-net) 
SNL is a SET-net together with a locality mapping @: Tr > N as in PTL-nets. 
Based on the semantics of the underlying SET-net and the semantics of nets 
with localities, one can then introduce four modes of execution, free-enabled 
/ min-enabled / max-enabled / lmax-enabled, in a straightforward way. 
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The modelling of a basic set membrane system as a SETL-net SNL(X) 
follows exactly the same lines as in Section 3 (note that all arcs have weight 
1 in this case, and there are no inhibitor nor activator arcs). Figure 7 shows 
the translation for the basic set membrane system ‘No in Figure 6. 


a Ti b c 
7 7% & ™ Ty 
Ox eee ia 
T12 
a 
c 
TT 
2 cae 
T21 a 
To 713 
aCe EO e@ @ 
122 5 


Figure 7: SETL-net SNL(o) modelling the basic set membrane system Yo. 


The tight correspondence between the membrane system % and the 
SETL-net SNL(X) is captured by the same two bijective mappings v and p, 
now restricted to sets, as in Section 3. Moreover, the key result (1) establish- 
ing the faithfulness of the translation obtained there can be re-stated for the 
current translation. This, in particular, means that the (finite and infinite) 
m-step sequences of SNL(*) faithfully represent m-computations of ©». 


5.2 Petri Net Analysis of Set Membrane Systems 


Moving from quantitative to qualitative membrane systems is an abstraction 
which may lead to a more tractable approach when it comes to answering 
vital questions concerning the evolution of systems. However, to take ad- 
vantage of this fact, the existing concrete analysis tools developed for the 
classical, quantitative, Petri net models need to be adapted for SET-nets. 
In [25], we have already made preliminary investigation into the synthe- 
sis problem which aims at an automatic construction of SET-nets exhibiting 
behaviour given in terms of a transition system. For set membrane systems 
this should contribute to insight in which evolution rules lead to certain 
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observed behaviour. 

By bringing qualitative (set rather than multiset) aspects to membrane 
systems, also interesting questions relating to expressive (generative) power 
emerge. For every mode, one can consider the possible evolutions of a system 
of aset membrane system (i.e., the computations of SETL-nets) as a language. 
These languages are regular subset languages. The study of subset languages 
of Petri nets was initiated in [40, 41] but still for the standard (quantitative) 
interpretation. There are a number of interesting theoretical questions and 
topics for the regular subset languages generated by SETL-nets under the 
four execution modes as well all regular subset languages. For example, one 
can consider: inclusion hierarchies; closure properties; and the complexity 
of equivalence/inclusion checking. Another group of problems here would 
be motivated by the target application area, i.e., biochemistry. For example, 
one can investigate: oscillatory behaviour (is it possible to have cycles from 
some point with at least /at most /specific evolution rules only); or vitality of 
the system (possible deadlock or partial death, i.e., some rules that can no 
longer be executed); or other state-related properties, like whether it would 
be possible for two different objects (types of molecules) to appear in a given 
compartment at some point together. 


6 Petri Nets and Extended Membrane Systems 


Basic (quantitative) membrane systems have over the past decade been ex- 
tended in several different directions, motivated either by their potential 
applications, or by their computational properties. For some of these ex- 
tensions, like catalysts and symport/antiport rules, there exist straightfor- 
ward translations to Petri nets (see, for example, |15]). For others, like i/o 
communication and rule creation/consumption, the correspondence between 
evolution rules and Petri net transitions is more involved, and the resulting 
nets are additionally equipped with inhibitor and/or activator arcs (see, e.g., 
[20]). 

First we consider an extended version of membrane systems for mod- 
elling how reactions may be triggered or blocked in the presence of certain 
molecules. The role of such molecules differs from that of catalysts which 
actively take part in reactions and are returned afterwards. An example is 
object a which acts as a catalyst in evolution rule 


r: {a,b} > {a, Dou} - 
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To model the subtle effect that the presence of molecules may have, mem- 
brane systems have evolution rules r of the form 


lhs” + rhs" | prom , inh” 


where pro” and inh” are multisets over V specifying respectively the pro- 
moters and inhibitors. The intuition behind pro” and inh” is that they test 
respectively for a minimal or maximal number of certain objects inside a 
compartment, but without consuming them. As a consequence, any number 
of rules can test for the presence of a single object at the same time. In 
order for r to occur there must be at least pro"™(a) copies of each symbol 
a in its associated compartment, and less than inh"(a) copies of each sym- 
bol a which occurs in inh”. Thus we retain all definitions introduced for 
basic membrane systems with only one change regarding the notion of a 
free-enabled vector multi-rule r. This is strengthened by additionally requir- 
ing that, for each 7 and r € r;, we have pro} < w; and, moreover, if a € inh” 
then w;(a) < inh"(a). 

PTL-nets are not expressive enough to model inhibitors and promoters 
because arcs between transitions and places indicate consumption and pro- 
duction of tokens (objects) rather than testing for their presence or absence. 
A possible way out is to use PTL-nets extended with range arcs [19]. Each 
such arc links a place to a transition and is specified by a closed interval (pos- 
sibly infinite) of non-negative integers. This interval indicates the range (a 
closed interval of natural numbers) for the number of tokens that should be 
present in the place to enable the occurrence of the transition. Clearly, like 
pro’ and inh”, range arcs can be used to model certain forbidden/required 
concentrations of molecules in a compartment. 

In [20], it has been shown that key properties of the modified transla- 
tion are very similar to those obtained in the basic case; in particular, the 
correspondence result (1) can simply be restated. Moreover, the treatment 
of causality developed for PTL-nets can also be extended (see [20]). 

What may come as a surprise, is that PTL-nets with inhibitor and ac- 
tivator arcs are also robust enough to model in a faithful way, membrane 
systems which have a dynamic structure due to rules which may thicken or 
dissolve membranes |21]. 

When it comes to the properties which might be expressed or investi- 
gated using reachability or coverability graphs, the situation changes dramat- 
ically when we move from PTL-nets to PTL-nets with range arcs. The reason 
is that the latter allow one to test for the absence of resources (zero-testing). 
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Nets with this kind of relationship between places and transitions (i.e., in- 
hibitor arcs) have been considered in [13]. The extension with inhibitor arcs 
gives the resulting model of Petri nets the expressive power of Turing ma- 
chines; net languages become recursively enumerable rather than recursive, 
and decidability for certain important behavioural properties, such as reach- 
ability and boundedness, is lost [13], partial solutions have been proposed by 
restricting the class of nets under consideration as, for example, in |4, 22]. 


7 Conclusion 


In this paper we described a number of results obtained while working at 
the interface between Petri nets (by now a classical formal model for deal- 
ing with distributed systems), and two recently proposed formal approaches 
aimed at dealing with computations inspired by biochemical reactions (i.e., 
membrane systems and reaction systems). Our overall experience was both 
illuminating and highly encouraging. We have found that the different mod- 
els share a number of important features which allowed us to compare them 
in a meaningful way (see, e.g., the faithfulness results (1) and (2)), but at the 
same time the original models did differ in some crucial aspects. The latter 
realisation has provided a direct motivation for importing concepts, such as 
the structuring of the molecules in a membrane system into compartments 
and the qualitative treatment of entities in reaction systems, into the realm 
of Petri nets. This has resulted in the introduction of new net classes, such 
as PTL-nets and SET-nets, which should be of a general interest to various 
application oriented research communities. For instance, nets with locali- 
ties together with the Imax-semantics are relevant for dealing with globally 
asynchronous locally synchronous (or GALS) systems. Moreover, we trans- 
ferred newly defined modelling concepts back into the domain of membrane 
systems, by developing set membrane systems. We therefore experience a 
harmonious synergy of the three, originally separate, models of computation. 

Results outlined in this paper open up a way to the adoption of Petri 
net analytical techniques in the areas of membrane computing and reaction 
systems. For example, one of the key advantages of Petri nets is that they 
support in a clear and unambiguous way concepts relating to causality and 
concurrency, and it turns out that the Petri net treatment of these concepts 
can be extended to membrane systems thanks to the so-called barb events 
which capture the intricacies of the lmax-semantics. We therefore feel that 
Petri nets are a robust model which can be suitably extended to provide 
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valuable insights into other, seemingly distant approaches. 
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